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ABELIAN HIGGS HAIR FOR BLACK HOLES 
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ABSTRACT 

We find evidence for the existence of solutions of the Einstein and Abelian Higgs 
field equations describing a black hole pierced by a Nielsen-Olesen vortex. This 
situation falls outside the scope of the usual no-hair arguments due to the non- 
trivial topology of the vortex configuration and the special properties of its energy- 
momentum tensor. By a combination of numerical and perturbative techniques we 
conclude that the black hole horizon has no difficulty in supporting the long range 
fields of the Nielsen Olesen string. Moreover, the effect of the vortex can in principle 
be measured from infinity, thus justifying its characterization as black hole "hair" . 
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Some years ago a collection of results were proved which established that the only long 
range information that a black hole could carry was its electromagnetic charge, mass and 
angular momentum. Thus, for example, lepton or baryon number were not good quantum 
numbers for black holes, despite being defined for a neutron star. This set of results came 
to be know as the "no-hair" theorems and, although not justified, gave rise to the much 
stronger, but for a while popularly held belief that the only non-trivial field configurations 
an event horizon could support were its massless spin one and spin two charges, Q,M, 
and J. Such a picture was not only misleading, but wrong; in spite of the structure of the 
original "no-hair" proofs, the no hair theorems only claim qualified uniqueness of static or 
stationary black hole spacetimes (see [1] for a discussion of no-hair folklore). 

The extrapolation of the no-hair folklore to include matter fields on the event horizon 
has been known to be false for some time, as has the extrapolation to include quantum 
effects. Black holes can be coloured [2-5], i.e. can support long range Yang-Mills hair. Such 
solutions are unstable [6, 7], thereby evading the usual no-hair uniqueness theorems, but 
they do nonetheless exist. Black holes also carry quantum hair [8, 9], which, although not 
locally observable, can be inferred via an Aharonov-Bohm interference in cosmic strings 
scattered on either side of the hole. Of course, the no-hair theorems are classical, but one 
might expect some modified argument to apply to quantum fields, so how is this long range 
hair mediated? More importantly, also at the semi-classical level, the hair has an effect on 
the thermodynamics of the black hole[10] caused by the phase shifts of virtual cosmic string 
loops dressing the euclidean horizon of the black hole. The existence of these euclidean 
vortices [11] seemed at first sight to be incompatible with some of the principles of the 
"no-hair" theorems, however, an examination of the appropriate theorem for the abelian 
Higgs model[12] revealed some assumptions which were not satisfied for the euclidean 
vortices[13]. Finally, it was shown that a small enough magnetically charged black hole 
would be unstable to the nucleation of an SU(2) t'Hooft Polyakov field configuration outside 
the horizon[14]. Technically of course, this is not new "hair" since the magnetic charge 
was already measurable at infinity, however, it illustrates nicely the difference between the 
"hair" and the short range massive fields which in this case can live outside and on the 
black hole horizon. 

It seems therefore, that one must be quite specific about what one means by hair on 
black holes, and we shall take the definition to mean a charge or property of the black 
hole measurable at infinity; whether or not fields can live on the horizon we shall refer to 
as dressing. The impression that a black hole horizon can support only a small number of 
long range massless gauge fields is somewhat misleading, an examination of the examples 
cited above indicates a common theme: when non-trivial topology is present in the field 
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theory, the situation is more subtle and dressing becomes a real possibility. In this paper 
we will examine the case of the abelian Higgs model, and show that the black hole can 
indeed sport long hair, namely a U(l) vortex. 

This paper was largely motivated by the work of Aryal, Ford and Vilenkin (AFV)[15], 
who wrote down a 'solution' for a cosmic string threading a black hole. More precisely, 
they wrote down an axisymmetric metric which consisted of a conical singularity centered 
on a Schwarzschild black hole. 

ds *= (l-^f) 1 dr 2 -r 2 de 2 -r 2 (1 - 4G^) 2 sin 2 Odcf (1.1) 

While such a metric was extremely suggestive that a true vortex would thread a black 
hole, it addressed only the gravitational aspect of the problem. In view of the delicacy of 
dressing the event horizon it seems necessary to analyze this part of the problem as well. 
The main potential obstruction to putting a vortex through a black hole is in having the 
vortex pierce the event horizon. Recall that the event horizon is generated by a congruence 
of null geodesies. If there exists a static vortex-black hole solution, then this congruence 
must remain convergence and shear free throughout the core of the vortex as it touches the 
black hole. This in turn translates to a relation on the stress-energy tensor, namely that 
Tq—T^ = 0. This is certainly true at the center of the string, where the energy and tension 
balance, but it is not clear that as we move away from the center and the null vector no 
longer aligns itself with the string worldsheet, that this balance will still be maintained. 
If a real vortex cannot puncture the event horizon of a black hole it raises questions as 
to how black holes and cosmic strings interact. Do they avoid each other entirely? Or 
does a black hole 'swallow up' a cosmic string caught in its gravitational grasp? And if a 
static equilibrium solution does exist, somehow avoiding the geodesic problem, how does 
the vortex pierce the horizon, and how does it circumvent the original no hair theorems 
for the Abelian Higgs model? 

In this paper we establish that the AFV solution can indeed be viewed as a thin string 
limit of some physical vortex by demonstrating that an Abelian Higgs vortex can thread the 
black hole. The layout of the paper is as follows: We first review the self-gravitating U(l) 
vortex in the next section. In section three, we examine the question of existence of the 
vortex in the black hole background, in the absence of gravitational back reaction. We find 
an analytic approximation to the solution for strings thin compared to the black hole, and 
present analytic and numerical results for strings of varying widths and winding numbers. 
We also examine global strings - a scenario for which we would not expect a gravitational 
solution - finding the reaction of the vortex to the event horizon to be different to that 
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of the local string. In section four we consider the gravitational back reaction of a single 
thin vortex, using the analytic approximation developed in section three. We derive the 
AFV metric, but find a subtle difference to their work involving a renormalization of the 
Schwarzschild mass parameter - all physically measurable results however agree. We also 
comment on the thermodynamics of the string-black hole system. Finally in section five 
we summarize and discuss our results, including an examination of more exotic systems 
such as a string terminating on a black hole - which has a gravitational counterpart. We 
also comment on the dynamical process of black hole cosmic string interaction and the 
second law of thermodynamics. Finally, we try to answer the question of whether the 
Abelian Higgs vortex is true hair, or whether it is just horizon dressing, and would be 
more appropriately dubbed a 'wig'. 

2. The Abelian Higgs Vortex. 

We start by briefly reviewing the U(l) vortex in order to establish notation and 
conventions. We take the abelian Higgs lagrangian 



where $ is a complex scalar field, -D M = + ieA^ is the usual gauge covariant derivative, 
and Fy, u the field strength associated with A^. We use units in which W = c = 1 and a 
mostly minus signature. For cosmic strings associated with galaxy formation rj ~ 10 15 GeV 
and A ~ 10" 12 . 

We shall choose to express the field content in a slightly different manner, and one in 
which the physical degrees of freedom are made more manifest. We define the (real) fields 
A, x and P^ by 



(2.1) 



$(x Q ) = rjX(x a )e ix ^ 
A^x a ) = -[P^x a )-V^ X (x' 



a 



)]■ 



(2.2a) 



(2.26) 



In terms of these new variables, the lagrangian and equations of motion become 



I A77 4 

C = ^V.IVX + V 2 X 2 P^ - j^F^F^ - ^-(A 2 - l) 2 



(2.3) 




(2.4a) 



V M F^ + 2e 2 rfX 2 P v = . 



(2.46) 
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Thus P M is the massive vector field in the broken symmetry phase, F^ v = V^P^ — V U P^ 
its field strength, and X the residual real scalar field with which it interacts, x is n °t in 
itself a physical quantity, however, it can contain physical information if it is non-single 
valued, in other words, if a vortex is present. In this case if the line integral of d\ around 
a closed loop (along which X = 1) is no n- vanishing, single- valuedness of $ then implies 

j> V^xcfa^ = [ x ] = 2tt?i for some n eH. (2.5) 

Continuity then demands (in the absence of non-trivial spatial topology) that X = at 
some point on any surface spanning the loop - this is the locus of the vortex. Thus the 
true physical content of this model is contained in the fields P^ and X plus boundary 
conditions on P M and X representing vortices. 

The simplest vortex solution is the Nielsen-Olesen vortex[16], an infinite, straight 
static solution with cylindrical symmetry. In this case, we can choose a gauge in which 

^ = V X (R)e^ ; P» = P Q {R)V ^ (2.6) 

where R = V\r]r c in cylindrical polar coordinates (r c ,<fi). The equations for X and P M 
greatly simplify to 

-Xd - § + ^ + \X,{Xl - 1) = (2.7a) 

-Po +^+/3" 1 X 2 Po = (2.76) 

where (3 = A/2e 2 = rn 2 scalar /ml ector is the Bogomolnyi parameter[17] {(3 = 1 corresponds 
to the vortex being supersymmetrizable) . Note that in these rescaled coordinates, the 
string has width of order unity. This string has winding number one; for winding number 
N, we replace x by Nx, and hence P by NP. 

This solution can be readily extended to include self-gravity by using Thome's cylin- 
drically symmetric coordinate system[18] 

ds 2 = e 2 ^\dt 2 - dr 2 ) - e 2 ^dz 2 - aV 2 ^ 2 (2.8) 

(where 7, ip, a are independent of z, (p) with the string energy-momentum tensor as the 
source: 

T^u = 2r, 2 V^XV u X + 2 V 2 X 2 P^P U - %F lur F v " - Cg^ (2.9) 
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in unrescaled coordinates. To rescale the coordinates, we set R = \f\r]r c , a = VXrja, 
and for future comparison, we write the rescaled version of the energy and stresses T ab = 
T ab /(\V 4 ) 

T ° = £ = e -2(7-^) X /2 + e * f + e 2( 7 -2^)p/2 /a 2 + ( ^2 _ 1)2/4 ^ ^ 

p a 2 

f| = -V R = -e-^-^X' 2 + € -4 t - e 2 ^- 2 ^P /2 /« 2 + (X 2 - l) 2 /4 (2.106) 
= e -2(7-V-) x /2 _ e —±-L- - e 2 ^- 2 ^P' 2 /a 2 + (X 2 - l) 2 /4 (2.10c) 



f z z = -V z = T °. 



Also for future reference, the Bianchi identity gives 

a' 

K + (v R - v+){- - v') + iv R + is = o. 

To zeroth order (flat space) 

a = R , ^ = 7 = 0, X = X , P = P , 

and (2.11) gives 

(i2P ofl )' = V 0(f> 

To first order in e = 8nGr] 2 the string metric is given by[19] 

-R 1 r R 



(2.10d) 
(2.11) 

(2.12) 
(2.13) 



a = 



l-e/ R(£ -V oR )dR R + e R 2 (£ - V oR )dR, (2.14a) 



7 = 2^ = e / RV oR dR. 
Jo 



(2.146) 



Where the subscript zero indicates evaluation in the flat space limit. Note that when the 
radial stresses do not vanish, there is a scaling between the time, z and radial coordinates 
for an observer at infinity and those for an observer sitting at the core of the string. The 
only case in which these stresses do vanish is when (3 = 1. In this case the field equations 
reduce to 

X' = XP/a 
P> = \a(X 2 - 1) 
a' = l-e[(X 2 -l)P + l] 

7 = ?/; = 



(2.15) 
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a first order set of coupled differential equations as one might expect from the fact that 
the solution is supersymmetrizable. 

We conclude this section by demonstrating the asymptotically conical nature of the 
corrected metric[19,20] (see also [21-24] for discussions involving model cores rather than 
vortices). Note that since the string functions X and P rapidly fall off to their vacuum val- 
ues outside the core, the integrals in (2.14) rapidly converge to their asymptotic, constant, 
values. Let 

t-R pR pR 

e / R{£ -V oR )dR = A, e / R 2 {S a -V oR )dR = B and e / RV oR = C (2.16) 
Jo Jo Jo 

then the asymptotic form of the metric is 

ds 2 = e c [dt 2 - dr 2 - dz 2 ] - r 2 (l -A + B lr c fe- c d(j) 2 

(2.17) 

= dt 2 - dr 2 c - dz 2 - f 2 (l - A) 2 e~ 2C d(j) 2 

where i = e°/ 2 t, z = e c l 2 z, and r c = e c / 2 (r c + B/(l — A)). This is seen to be conical 
with a deficit angle 

A = 2tv(A + C) = 2yre J R£ dR = 16w 2 G J r c T°dr c = 8ttG^ (2.18) 

where \i is the energy per unit length of the string. Notice that the deficit angle is inde- 
pendent of the radial stresses, but that there is a red/blue-shift of time between infinity 
and the core of the string if they do not vanish. 

3. String in background Schwarzschild metric. 

In solving on a Schwarzschild background there are two coordinate systems we could 
consider. 

i) Spherical (Schwarzschild) coordinates 

ds 2 =(l- dt 2 -(l- 1 dr 2 s - r 2 (d6 2 + sin 2 Odcf 2 ) (3.1) 

This is the usual Schwarzschild metric. This coordinate system is good for analyzing the 
existence of a background solution, but, not being tailored to the symmetries of the full 
problem, it does not deal with gravitational back-reaction well. 

ii) Axisymmetric (Weyl) coordinates [25] 

Here 

^2 _ fli+fla -1GM ul (Ri+R 2 +2GM) 2 / ,2 , j2\ 2 (fli+fl 2 +2GM) ,,2 

ab ~ R 1 +R 2 +2GM Ui 4R 1 R 2 c) ' c (R 1 +R 2 -2GM) U ^ y ^' 
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where 

R\ = (z - GM) 2 + r 2 , R\ = (z + GM) 2 + r 2 c (3.3) 
The transformation between the two systems is given by 

z=(r a - GM) cos 9 , r 2 c = r s (r a - 2GM) sin 2 9 (3.4) 

Weyl coordinates are appropriate for analyzing the gravitational back reaction, however 
are rather cumbersome for the background solution problem. 

The rest of this section is devoted to arguing existence of a vortex solution in a black 
hole background. Formally, this means taking the somewhat artificial limit Grj 2 — > 
keeping GM fixed. It is straightforward to show that it is consistent to take 

P„ = PV^ (3.5) 

and letting 



r = V\rjr 3 
E = V\ V GM 

in the Schwarzschild metric, the equations of motion for X and P are 



(3.6) 



-la r [r(r - 2E)d r X] - -^^[sin^X] + \X{X 2 - 1) + = (3.7a) 



sin 9 



d r [(l - 2E/r)d r P] + — r d e {Qsc9d e P] - /3~ 1 X 2 P = (3.76) 

We proceed at first by taking a "thin string limit", in other words we assume 2E ^> 1. 
This is equivalent to requiring M ^> 1000 Kg for the parameters of the GUT string. We 
will then consider thicker and higher winding number strings. 
First of all, let us try 

X = X(rsm9),P = P(rsin9) (3.8) 

substituting these forms into (3.7), and writing (suggestively) rsin6* = R we get 

-X" - X'/R + \X(X 2 - 1) + XP 2 /R 2 + 2E sin 2 9/r[X" + X'/R] = (3.9a) 
P" - P'/R - P~ 1 X 2 P + 2E sin 2 9/r[-P" + P'/R] = (3.96) 

where prime denotes derivative with respect to R. If we were to input X, P = Xq, Pq then 
these equations would be satisfied up to errors of the form 

2£sin 2 # r , . . . /0 
x [other terms in equation] (3.10) 



However, since rsin6> = R ~ 0(1) in the core of the string, sin6> = 0(l/r), and the 
errors are 0(E/r 3 ) < 0{1/E 2 ) < 1. Thus the Nielsen-Olesen solution is in fact a good 
solution throughout and beyond the core of the string, whether near the event horizon 
or not. By the time the premultiplying term in the errors is significant, we are well into 
the exponential fall-off of the vortex, and essentially in vacuum. However, since we are 
interested in showing that the event horizon can support the vortex, for completeness we 
include the solution to order 1/E on the horizon: 

l-X~e- 2Ee + e - 2E ^- \ P ~ e-WF 7 * + e -2E^(n-9) (3n) 

Of course, these are only approximate forms, and do not prove existence of a solution to 
(3.7). However, they will provide a good approximation to the true solution, if such can 
be shown to exist. We will provide such evidence in the form of numerical solutions later 
in this section. For the moment, we conclude this description of the thin string limit by 
examining the equations of motion in the extended Schwarzschild spacetime in terms of 
Kruskal coordinates. 

One problem with using Schwarzschild coordinates for our analysis is that they are 
singular on the event horizon. This is, of course, purely a coordinate singularity, but since 
demonstrating a dressing of the event horizon is central to this paper, we will examine 
the thin vortex solution in coordinates that are not singular at the event horizon, namely 
Kruskal coordinates, to convince the reader that the analytic approximation really does 
hold true at the event horizon. Kruskal coordinates are based on the incoming and outgoing 
radial null congruences of the Schwarzschild spacetime, but we shall instead use a Kruskal 
'time' and 'space' coordinate defined as follows. 



T=(r s - 2GM) 1 /V*/ 4GM sinh [j^j 
Y = (r s - 2GM) 1 / 2 e r °' AGM cosh ^ 



4GM 



(3.12) 



In terms of which the metric is 



ds 2 = 16G2M \ -^/ 2GM (dT 2 - dY 2 ) - r 2 (d6 2 + sin 2 Odf) (3.13) 

This metric is clearly regular away from r s = 0, and the future and past event horizons 
are represented by T = Y > and — T = Y > respectively. What we would like to 
show is that X = X (rsm6) and P = P (r sin 9) are good solutions (expressed in Kruskal 
coordinates) in the vicinity of the horizon, \T\ = \Y\. 
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First note that near the horizon 



r = VXqr s = 2E - e^y/X^T 2 - Y 2 ) + 0(T 2 - Y 2 f 



(3.14) 



where here e = 2.718... is the natural number. This implies that 



(3.15) 



and hence 



-^d T g TT ^d T + -^d R g RR ^d R 



X = Xrf 



{2E r Km 2 OX'\R)+ { -^—p-X'(R)sme 

2E 



2E 



i a g TT V9 8 i 8 g RR V9 8 ■ 

~^ d T 9 . 2 n OT + —OR 9 . 2 n dR 

y/g rj sin 9 y/g rj sin 9 



P = \rf 



(2E-r) p// _ P' 



2E 3 



E 2 R 



(3.16) 



In other words, as in Schwarzschild coordinates, the Nielsen- Olesen solution solves the 
equations of motion to 0(E~ 2 ) near, on, and even beyond the event horizon. Indeed, 
replacing 

r = v / A?]r*" 1 {2GMlog(y 2 -T 2 )}, (3.17) 

where r*(r s ) is the tortoise coordinate, indicates that the approximation holds true well 
within the event horizon for larger black holes. 

Having established that it is possible for the horizon to support N = 1 vortices, we 
now turn to the large iV case, where analytic approximations are available [26]. First of 
all, note that a string with N ^> 1 has a core radius of order y/~N. To see this, consider 
the rather special case of j3 = 1 (for general (3, see [26]). Setting X = £ N , we have from 
(2.15) 

-2N _ -J ^ _ ^_ 

t R 



np' = \R{e N - 1) 



(3.18) 



in the absence of gravity. These are in fact the same as the large-N equations for general 
/?, it is in the sub-leading terms that the two cases differ. These are solved in the core by 



P = 1 



4N 



■ £ocRe- R2 / 2N 



(3.19) 



using < £ < 1. The transition to vacuum, and hence a different approximate solution to 
the above can be seen to occur quite abruptly at R = 0(y/~N). Thus the condition that 
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the string be thin compared with the black hole is now 2E 3> y/N, and in that case the 
previous arguments still apply, since 



2E sin 2 6 _ 2ER 2 2EN N 

" - -73- < -73- ^ < ' 

i.e., the Nielsen Olesen solution for a large-N string is good also on the event horizon. 

Now consider the opposite limit, where the string is much bigger than the black hole, 
VN 3> 2E . The black hole sits well inside the core of the string, in a region where 
I K 0, P 1 — pR 2 . In the absence of a black hole the P-equation simply states that 
the magnetic field is constant throughout the core of the string. Notice that since we are 
ignoring X 2 in the large-N expansion, (3.9b) implies that the presence of the black hole 
does not affect the large- N P-equation, so we still find P(R) ~ 1 — pR 2 . The magnetic 
field will still be constant (and equal to — 2p) in the string core. However, its value may 
change due to the black hole. Notice that we may expect a slight "squeezing" of the string 
core due to the black hole. To see this, consider rewriting the P-equation (3.9b) in the 
form 

P" - (P'/R) - /T 1 (r, 9)PX 2 = 0. (3.20) 

On the equatorial plane of the black hole, sin# = 1, and (3~ 1 (R) = /3 _1 /(1 — 2E/R). For 
R 3> 2E we have the vacuum solution X = 1, P = 0. As we come in towards the horizon, 
P has to leave its vacuum value - however, the effective value of (which measures an 
"effective mass" for P) is increasing. Compared with the situation where there is no black 
hole, P should be more reluctant to leave its vacuum value. The magnetic field will remain 
zero for as long as possible and, as a result, the string core is somewhat smaller around 
the black hole. Note that this argument does not apply for global strings, where numerical 
simulations indeed show that it is a much smaller effect. 

Now consider the X equation, to leading order in N (r > 2E): 

i -^)(?H-W^)(?) , -£ +o « <-> 

We want to show that the solutions to this equation are regular at the horizon. The 
equation becomes singular at sin6> = 1, r = 2E, (or z = 0, R = 2P), the equatorial plane 
of the horizon, so let us integrate the equations to leading order in 1/N and z/R: 



(3.22) 
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We find 



.l(i-|p£ 2 ) 



e - P [^]y/R(R-2E) ^ 



i = K R-E + \/R(R - 2E) 



(3.23) 



which is finite at R=2E. The constant K can be fixed by the requirement that £ ~ 1 when 
P ~ 0, i.e., at R = 1/y/p. 

Finally, note that the horizon seems to be capable of supporting global strings as well, 
in spite of the fact that we do not expect to be able to account for gravitational back 
reaction consistently, since the energy per unit length of a global string diverges and the 
self-gravitating global string spacetime is singular [27, 28]. To find the global string solution 
in the black hole background we simply set P = 1 everywhere, to find 



-X" - X'/R + \X(X 2 - 1) + NX/R 2 + 2£sin 2 9/r[X" + X'/R] = (3.24) 



The case where the global string is thin compared to the black hole works as before - the 
vortex is essentially undisturbed. In the case where the string is bigger that the black hole 
we can take p — > in (3.23), so the solution is again regular at the horizon. We conclude 
that the presence of the black hole is, if anything, less noticeable than in the local string 
case, as can be seen in the numerical simulations described next. 

Numerical results 

We will now provide confirmation of the previous analytic arguments by means of a 
numerical solution of the equations of motion outside the event horizon. To this end, we 
note that the equations are elliptic outside the event horizon, parabolic on it and hyperbolic 
inside it. Some care is therefore required with specification of the boundary conditions. 

At large radii we want to recover the NO solutions, while the symmetry axes outside 
the horizon must form the core of the string: 




(3.25) 



On the horizon the equation turns parabolic, taking the form 




r=2E 



AE 2 sin 9 



1 



d e [sm0d e X} + \X{X 2 - 1) + 



4£ 2 sin 2 9 ' 



N 2 XP 2 



(3.26a) 




sin6> 



de[csced 9 P] -I3~ 1 X 2 P. 



(3.266) 



r=2E 



AE 2 
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A given choice of X(2E, 9) and P(2F, 9) on the horizon therefore directly implies d r X{2E, 9) 
and d r P(2E, 9). Generally boundary values and normal derivatives cannot be specified si- 
multaneously when solving an elliptic equation: either suffices to determine the solution 
completely. A particular choice of X and P on the horizon and the corresponding normal 
derivatives will not normally generate the same solution outside the horizon. The bound- 
ary conditions on the event horizon are thus determined by the requirement that both 
solutions coincide. 

A practical algorithm for solving the equations of motion in a Schwarzschild back- 
ground numerically then is as follows. We employ a uniformly spaced polar grid {(r^, 9j)}, 
with boundaries at r = 2E and a large radius r/, > 2E, and 9 ranging from to ir. Then 
we approximate the derivatives with finite-difference expressions on the grid. Writing Poo 
for the value of the field F at the grid point (r^, 6j), and similarly F± for F(r i ± 1 , 9j) and 
F ± for F(ri, 0j±i), we obtain the finite difference equations 

2 (-[ _ E\ X +0 -X_o cot e X 0+ -X _ (-. _ 2E_\ X +0 +X_ X Q+ +X ^ 
X — r \ L r) 2Ar r^ 2A9 t I 1 r ) Ar 2 [ r 2 Aff 2 /q 27 a ) 

°° h 2E \ 2 I 2 i l(v2 -in | ( NP 00 \ 2 ' ' 

I 1 ~~ "J "Ar 1 + ^AP" 2~l A 00 - J-J + 



P+o-P-o , nPo±-Pq-_ , A 2^ P+q+P-q , P 0+ +Po- 
_ cot y_ sg _ 3 _ + ^ _ _j h r2Ag2 



1 2Ar 

00 = 



2 

00 



(3.276) 



inside the grid, and 



Xoo = - Ar , 2 ± e2 __„ - 4A J„ 72 i (3.27c) 



E 
Ar 



+ Ak+^0 2 -l) + |(S)' 
P 0+ +Po- _ t g P 0+ -P - 



(3.27d) 



on the horizon. 

Initial values for X and P are assigned on the boundaries according to (3.25); on the 
horizon we initially set X = 0, P = 1. X and P are then iteratively adjusted on the 
interior grid points according to (3.27a,b), analogous to the Gauss-Seidel scheme for linear 
elliptic equations [29]. After each pass through the interior grid points, the r-gradients of 
X and P just outside the horizon are calculated, and equations (3.27c,d) iterated to derive 
new values for X and P on the horizon (for given r-gradients, the equations on the horizon 
are one-dimensional elliptic equations) . The whole process is then iterated to convergence. 
In order to speed up convergence, the grid is over-relaxed: instead of replacing X and P 
by the right-hand sides of eqs. (3.27), wX new + (1 — w)X with 1 < w < 2 is used. The 
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optimal value for the over-relaxation parameter w is found by trial and error, and depends 
on the number of grid points and on the differential equation. 

Sample results are presented in Figures 1-7, and confirm the analytic arguments above. 
Figures 1-4 show a sequence of solutions with increasing winding number (and therefore 
string thickness) threading an E=10 black hole. Qualitatively, the string simply continues 
regardless of the black hole, though some mild pinching of the magnetic flux does take 
place. Figures 5 and 6 compare a local and global string with the same winding number 
and 'width'; the global string is apparently fatter due to the power law, as opposed to 
exponential, fall-off in the fields. Figure 7 shows a comparison between the numerically 
obtained solutions and the Nielsen-Olesen analytic approximation that will be used in the 
next section. As can be seen, even for a string of non-negligible thickness, this is still an 
excellent approximation. 

4. Gravitating strings 

In order to get the gravitational effect of the string superimposed on the black hole, 
we need to consider a general static axially symmetric metric 

ds 2 = e 2 ^dt 2 - e 2 ^ (dz 2 + dr 2 c ) - a 2 e"^# 2 (4.1) 

where ^,7,0; are independent of t, 0. Notice that this is related to (2.8) through z — > 
it, t — > iz. We then apply an iterative procedure to solving the field equations, starting 
with the background solutions (3.2) and the Nielsen-Olesen forms of X and P, and ex- 
panding the equations of motion in terms of e = 8nGr] 2 which is assumed small, (e < 10 -6 
for GUT strings.) We first rescale coordinates to bring them into line with the rescaled 
Schwarzschild coordinates used in the previous section 

p = V\rjr c 

C = V\ V z (4.2) 
a = v^Xrja 

and rewrite 

Rl = ((-E) 2 + p 2 

(4-3) 

R 2 = (( + E) 2 + p 2 

In terms of the rescaled coordinates and energy-momentum tensor, the Einstein equations 
become 

a a + a pp = -ev^(f c c + Tfl (4.4a) 
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(«V>,c),c + W,p),p = ^v^Cft - j< - f; - tJ) 

(« 2 C + " 2 p )7,p = ey/=g(a jP f£ - a C T C P ) + aa iP (^ - V> 2 C ) 
+ 2aa )C -0,cV',p + a ,p a , PP + a ,C a ,Cp 

(« 2 C + « 2 pHc = -ev^Ca.C^ - ct, p f£) - aa,((^ p - Vj) 
+ 2aa^^, P ~ a ,C a ,pp + a ,p a ,(p 
%pp + 7,CC = ~ V? C ~ ee^-^ft 

where the energy momentum tensor is given by 



(4.46) 
(4.4c) 

(4.4d) 
(4.4e) 



, x ^ 2 ^ 2 
Y 2 P 2 



f c c = V(X) + 



+ 



rpp _ _2g-2(7— 0) 



+ (* , + * c) 



X 2 P 2 

a 2 e -2i, 

X 2 P 2 

a 2 e~ 2 ^ 



1 ' "' + {X 2 - XI) 



(P 2 - P 2 



(P 2 c - P 2 p) 

P~l a 2 e -2^ 
P tP 

fi^^e- 2 ^ 



-2( 7 -V) 



-2( 7 -V) 



-2( 7 -V) 



,-2( 7 -V) 



(4.5) 



We now write a = ao + ecti etc., and solve the Einstein equations (4.4), and the string 
equations (2.4), iteratively. To zeroth order, we have the background solutions 



«0 = P , V>0 = \ log l\X^-+2%M . 70 

X = X ( J R) , P = Po(R) , 



= ^log 



(R 1 +R 2 -2GM)(R 1 +R 2 +2GM) 
4R 1 R 2 



(4.6a) 
(4.66) 



where i?i and i? 2 were defined in (3.3). In these coordinates = rsin6> = pe~^°, so (4.6b) 
indicates that many of the terms in T 6 a are simply functions of R. 

Before proceeding to calculate the back reaction however, it is prudent to check that 
the energy-momentum tensor (4.5) will admit a geodesic shear free event horizon. Recall 
that we require Tq — T£ = on the horizon in Schwarzschild coordinates. This is clearly 
satisfied at 9 = 0, where the energy and tension balance, but what about 6> 7^ 0? In Weyl 
coordinates, this corresponds to Tq — T? = for p — > 0, C From (4.5) we see that 
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this is given by 



(3P 2 



ae -2i, ^ ,P 



8R\R 2 



(#1 + R 2 + 2EY \Jp 



(3P'(R ) 2 
R 2 



(4.7) 



+ X'(Rf 



All terms in this expression remain finite and non-zero as p — > (R — > 0) except for 
Using the transformation (3.4) between the Schwarzschild and Weyl coordinates, we have 



OR _ p{r-E sin 2 9) 
dp R ± R 2 sin 9 
dR Er sin 9 cos 9 

~ hji 2 



(4.8) 



hence ^ -> as p -> and f ° - fp is indeed zero on the horizon. Thus there is no 
gravitational obstruction, at least in this linearized method, of painting the vortex onto 
the horizon. 

Proceeding with the analysis of the energy momentum tensor, we use R\ — R\ = AEC,, 
which gives £1^2 = E 2 sin 2 9 + r(r- 2E), thus near the core of the string, where sin 9 = 
OiE- 1 ) 



R 2 + ri = (i- 2 JL sin 2 e] ~ = e 2 (^o). 

,P R X R 2 \ r J R X R 2 



(4.9) 



Therefore, in and near the core of the string the zeroth order rescaled energy-momentum 
tensor now reads 



Y 2 P 2 

TS, = n*o)+ ° 



+ 



p/2 



+ X' 2 + 0(E~' 2 ) 



R 2 I3~ 1 R 2 



v2 p2 p/2 

V(X ) - - + X' 2 + 0(E- 2 ) 



R 2 P^R 2 



X 2 P 2 



K = v(x ) + ^ + 



R 2 



i* = V(X ) + 



R 2 



p/2 

ii-^R 2 

p/2 



12 



+ *0 



12 



(i-^R 2 



f o P c = -2e~ 2 ^o-i>o)R c R 



p/2 

+ A o 



fi-^R 2 
16 



{R\-R\)e~^-^ 
+ 0(E~ 2 ) 



(4.10) 



and the combinations used in equations (4.4a,b,e) are all purely functions of R. This 
strongly suggests looking for metric perturbations as functions of R. However, we must 
check that the left hand sides of these equations can be written as appropriate functions 
of R. 

Consider R = pe~^° , 



R 



R 



and 



imply 



R, P = {- p -^ Q , P )R 



R PP = -(f"/> ,p),P^ - R A- n - %) 

P p K 



(4.11) 



(4.12) 



^,CC + R ,pp - 



R 2 C + R 2 P R 



R 



R 



Esin9 



rR x R 2 

where we have used the zeroth order equation of motion for ipo. Therefore 

d 2 



0(E- 3 )e 2 ^-^ o ) 
(4.13) 



d 2 + d 2 p = (R 2 C + Rl)^ 2 + (R XC + Rm ,)± = ftoo^o) 



dR 



dR 2 



+ 0{E~ 6 ) 



(4.14) 



in the core of the string. Exterior to the core, the vacuum equations will apply. We now 
solve (4.4a,b,e) to first order in e, namely 



<*i,CC + a hPP 



-2pe 2 ^°\V(X ) + ^-) 



R 2 



= -pe 2 <™-+ o) (£ -P 0R ) 



ai,C0o,c + "l.pV'o.p + P^hCC + (pV>i,p),p = pe 2(7 ° * o) 



p/2 

%-vix Q) 



7i,« + 7i,pp + 2Vo>i, P + 2Vo, C ^i,c = -e 2 ^-^f* 

— g2(7o— 0o)p 



(4.15a) 



(4.156) 



(4.15c) 



where £ and the Vs are given by (2.10). 

We first solve for a±. Note that there is an a = p in the 
This suggests that we write 

cki = pa(R). 



^on the RHS of (4.15a). 



(4.16) 
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a(R) then satisfies 



a"(R) + ^a'(R) = -[E -V 0R ] (4.17) 



Thus 



< R ) = - j jp J R 2 [^-v 0R ]dR 

= - J R[£ - V 0R ]dR + i J R 2 [S - V 0R ]dR (4.18) 

this is readily seen to have the asymptotic form 

a{R)~-- + ^- (4.19) 
e eR 

(where A,B are given by (2.16)) and solves the vacuum equations. 

Setting ipi = tjji(R), and using the form of ol\ given by (4.16,18), we see that (4.15b) 
becomes 

V>" (R) + ^i(R) = hPoR + ( 4 - 20 ) 

which is solved by 

Vi = \ f\ f R(V oR + V o4> ) = \ J RV oR (4.21) 

using the zeroth order equations of motion (2.12). Thus ipi tends to a constant (C/2e) at 
infinity which is also a vacuum solution. 

Finally, setting 71 = 71 (R), and using the form of ipi given above, (4.15c) reduces to 

li(R) = Vo* (4-22) 

and 

7i - || V 04> dR = J RV 0R dR = 2^i. (4.23) 

Thus, the corrections to the metric written in this form are almost identical to the self- 
gravitating vortex solution. In fact, using these corrections, we see that the asymptotic 
form of the metric given by 

ds i _ e c [e 2 ^dt 2 - e 2 ^-^\dr 2 c + dz 2 )] - r c 2 (l - A + ) 2 e" c e" 2 ^# 2 (4.24) 
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in the Weyl metric, or 



eC / _ 2GM\ f2 _( 1 _ 2GM\ 1 dr 2 r 2 d6 2 ] _ r 2 {1 _ A+ I3 fc-c^^ 

V r s J V r a J V\ V r s sme 

(4.25) 

in the Schwarzschild metric. Note that although the B-term appears to distort the event 
horizon, B/y/\r]r s = 0(Gfi) x 0(E~ 1 ) and hence represents an effect outside the regime 
of applicability of our approximation. We therefore drop this term, rescale the metric so 
that time asymptotically approaches proper time at infinity, t = e c l 2 t etc. and setting 
M = e c l 2 M we have 

ds 2 =(l- dt 2 -(l- dr 2 - r 2 d9 2 - r 2 (l - A) 2 e~ 2C sin 2 9dc/> 2 (4.26) 

We thus see that our spacetime is asymptotically locally flat with deficit angle 2n(A + 
C) = 8nGfi. Thus by using a physical vortex model, we have confirmed the results of 
AFV. However, note that the presence of the radial pressure term e~ c has modified the 
Schwarzschild mass parameter at infinity to M = e c l 2 M. The gravitational mass of the 
black hole has therefore shifted to 

M g = M = e c/2 M (4.27) 

The inertial mass of the black hole, or its internal energy, can be found by considering the 
black hole as being formed by a spherical shell of matter infalling from infinity. Due to the 
deficit angle this has mass 

Mi = M(l — A)e~ c = M g (l — 4G/i) (4.28) 

thus the inertial mass of the black hole is actually less that its gravitational mass, however, 
since we cannot accelerate the black hole without accelerating the string, it is perhaps more 
correct to refer to this as the internal energy of the black hole. We conclude this section 
on the gravitating string-black hole system by remarking on the thermodynamics of the 
system. 

Either by euclideanization, or by considering the wave function of a quantum field 
propagating on the black hole background, one can see that the temperature of the black 
hole is T = (3~ x = l/8nGM g . We denote the thermodynamic quantity T _1 as f3 to 
distinguish it from the Bogomolnyi parameter, additionally, we have set the Boltzmann 
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constant, k, to unity. Since the spacetime is no longer asymptotically flat, euclidean argu- 
ments must be interpreted with care; nonetheless, by a somewhat non-rigorous partition 
function calculation we confirm the AFV result that the entropy of the string-black hole 
system is 

s = ^ 1 - A ^° = w A (4 ' 29) 

Thus although the temperature of the black holes is unchanged in terms of the gravitational 
mass measured at infinity, and although the area-entropy relationship is unchanged, since 
the internal and gravitational masses are no longer equal, the entropy of the black hole 
with the string is less than that of a black hole of the same temperature (i.e. gravitational 
mass) without the string. 

It is interesting to use these thermo dynamical results to examine the dynamical situa- 
tion of a cosmic string - black hole merger. If one demands that the gravitational mass of 
the black hole is fixed, then the temperature of the black hole remains unchanged, but its 
entropy decreases. If one demands conservation of internal energy, then the temperature 
decreases and the entropy increases. Clearly, thermodynamics indicates that conservation 
of internal energy is the correct condition to use. It is interesting to note that in this case 
the change in gravitational mass is 

5M g = Mg — Mi = 2 x 2GM g x fi (4.30) 

i.e. the change in gravitational mass is equivalent to the length of string swallowed up by 
the black hole as seen from infinity times its energy per unit length. In this sense, the 
vortex at infinity is direct hair, conveying exact information as to the last 4GM/I units 
of matter that the black hole swallowed. We will take up this theme further in the next 
section. 

5. Summary and Discussion 

In this paper we have provided evidence, both analytical and numerical, that U(l) 
abelian Higgs vortices can pierce a black hole horizon. We have shown that there is 
no gravitational obstruction to this solution, and at this point it is perhaps worthwhile 
detailing how our solution avoids the revamped abelian Higgs no-hair theorem[30]. A 
simple answer would be that the string system is not spherically symmetric, however, many 
of the steps in [30] can be generalized to include more generic situations. Indeed, recent 
interesting results of Ridgway and Weinberg[31], who show non-spherically symmetric 
dressing (although not hair) of black hole event horizons, indicates that spherical symmetry 
should not be a prerequisite of no-hair theorems. The main reason our solution evades 
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such a no-hair 'proof is that a vortex mandates a non-zero spatial gauge field, which then 
destroys the inequalities on which the no-hair results are based. In particular, the fact that 
the field P^ has lines of singularity (corresponding to the vortex cores) explicitly breaks 
the argument given in [30] as to the vanishing of Pj. 

But is the vortex dressing or hair? The thermodynamical argument seems to indicate 
that it is hair - telling us about the last AGMji units of mass the black hole swallowed. 
Can we make this argument stronger? Suppose instead of considering a vortex threading a 
black hole we consider a single vortex terminating on a black hole. This has a gravitational 
counterpart in the guise of a uniformly accelerating black hole connected to infinity by a 
conical singularity [32], therefore we can ask whether there exists a particle physics vortex 
counterpart to this setup. One immediate difference with the previous situation is that the 
metric here is non-static, however, that can be remedied by introducing a second black hole 
attached to infinity by a second string placed so that its gravitational attraction neutralizes 
the uniform acceleration[33]. This leaves us with the topological question of how to paint a 
single semi-infinite vortex onto a black hole event horizon. Recall in section two, when the 
transformation to the real variables X and P^ was performed that the phase of the Higgs 
field, x, was purely gauge, and only acquired physical significance via boundary conditions 
on Pp. It was remarked that if spatial topology was trivial, any surface spanning a loop 
with a winding of x would have a vortex on it. If the space is non-trivial, then the question 
of whether any surface spanning such a loop must have a vortex reduces to a question of 
topology: if the first Chern class of the U(l) bundle is trivial, then any spanning surface 
must have a vortex. 

Recall that in Kruskal coordinates, the extended Schwarzschild spacetime contains a 
wormhole: the t =constant surface. This has topology S 2 xIR with two asymptotically 
flat regions. Thus the spatial topology of the Schwarzschild black hole is non-trivial, and 
is homotopically equivalent to a sphere. The issue of whether a vortex can terminate on 
a black hole therefore reduces to that of placing vortices on 2-spheres - a well studied 
problem (see [34] and references therein). In our case, the answer is to take two gauge 
patches, e.g. 

h = {e,(t)\e <tt/2 + s} 
h = {e,(f)\e>7r/2-s} 

for any 5 < tt/2. Then define the (gauge) transition function 

912 = e* 
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such that 

$1 = #i2$ 2 ; A llt = A 2m - ^-d^ 12 

on the overlap, and we take $ 2 = ^2 = 0. Thus we have a vacuum on the southern 
hemisphere and a vortex on the northern hemisphere connected via a non-singular gauge 
transformation on the overlap. Of course, if this two sphere could be shrunk to zero radius, 
then this would not be an allowed gauge transformation, but since the two-spheres in the 
Schwarzschild spacetime have a minimum radius, 2GM, there is no topological obstruction 
to this definition, and we can therefore have just a single vortex connected to the black 
hole. In terms of the extended Schwarzschild spacetime, this vortex enters the black hole 
via the North Pole, goes down the wormhole, and emerges from the North Pole of the 
black hole in the other asymptotic regime. The string worldsheet itself looks like a two- 
dimensional black hole, but occupies only the 6 = portion of the full four-dimensional 
Penrose diagram. We have not verified that the Nielsen- Olesen solution can be painted on 
to the non-static accelerating black hole spacetime, however, based on the static evidence 
and the lack of a topological obstruction, it would be very surprising if it could not be. 

This now leaves us with the question of how a black hole might have got just a 
single semi-infinite vortex in the first place. Certainly it cannot happen as the result of 
interaction between an infinite vortex and a black hole, so let us consider what the presence 
of the vortex actually means. When a single vortex is present on the two-sphere, more 
than one gauge patch is necessary for a non-singular description of the physics. This is 
analogous to the Wu-Yang[35] description of a Dirac monopole. Indeed, given that in each 
case we are dealing with the same mathematical object (a U(l) bundle over the sphere) 
the only real difference between the two cases is the spontaneously broken symmetry. 
Thus the interpretation of the vortex is that it is localized 'magnetic' flux emanating 
from the black hole. In terms of the dynamics of phase transitions in the early universe, 
one is led to a picture of a magnetically charged Reissner-Nordstrom black hole prior 
to the phase transition having its flux localized in the vortex after the phase transition. 
Thus the information (namely 'magnetic' charge) which one would not normally expect 
to be able to measure corresponding as it does to a massive field, is indeed preserved for 
external observers to see in the form of the long vortex hair stretching to infinity. We can 
correspondingly imagine a charge 2 Reissner-Nordstrom hole becoming a Schwarzschild 
hole with two vortices extending to infinity from its opposite poles, which would then be 
of the AFV form described in section four, where it was the energy-momentum rather than 
the orientation of the vortices that was relevant. 

In the light of this evidence, we claim that the Abelian Higgs vortex is not simply 
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dressing of the black hole, as the SU(2) monopole is, but is true hair - carrying information 
from the black hole to infinity. 

Note added in proof. After this work was completed, we were informed that Eardley 
et.al. [36] had also developed the gauge patch description (Sec. V) in order to argue the 
instability of NO vortices to black hole nucleation. Additionally, the conjecture in Sec. V 
that the NO vortex could be painted on to the nonstatic C-metric has since been verified 
in [37]. 
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Figure 1. Numerical solution of the Nielsen-Olesen equations with iV = l,/3 = l/2ina 
Schwarzschild metric (E = 10) background. The event horizon is indicated by a semicircle. 
Evidently the presence of the black hole horizon hardly affects the string structure at all. 
This solution was calculated with 100 radial and 100 azimuthal grid points, out to radius 
r L = 60. 
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Figure 2. As figure 1, but for winding number 5. 
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Figure 3. As figure 1, but for winding number fOO. The string is noticeably pinched. 
(tl = 100 for this calculation.) The undulations in the outer contours occur on the scale 
of one grid cell, and are an artifact of the contouring package's conversion from polar to 
rectangular coordinates. 



~ i 1 1 r 



b 

II 

PL, 



TO 
Sh 

o 

a 
o 
o 

I 

PL, 



05 



n 1 1 r 





z 



0Q- 




27 



Figure 4. As figure 1, but for winding number 400. The event horizon is now entirely 
inside the core of the string, which is slightly pinched, (rj, = 150) 
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Figure 5. As figure 1, but with E = 1, N = 1, tl = 15. String and black hole have 
comparable radii, but distortion of the string by the background is still rather mild. 
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Figure 6. As figure 5, but for a global string. In agreement with our analytical arguments, 
the effect of the black hole on a global string is weaker than on the corresponding local 
string. 
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Figure 7. Illustration of the relatively small effect of a black hole horizon (E = 10) on 
a local (j3 = 0.5, N = 50) string. In these panels, solid lines show values of the field in 
the black hole background, and dashed lines the values at corresponding positions in a flat 
metric. Upper and lower panels show X and P, respectively, while left and right panels 
show cuts along the equator (9 = rc/2) and around the horizon (r = 2E). 
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